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Abstract 



We consider the (smoothed) average correlation between the den- 
I sity of energy levels of a disordered system, in which the Hamiltonian 

is equal to the sum of a deterministic Hq and of a random potential if. 
Remarkably, this correlation function may be explicitly determined in 
the limit of large matrices, for any unperturbed Hq and for a class 
of probability distribution P{ip) of the random potential. We find a 
compact representation of the correlation function. From this repre- 
sentation one obtains readily the short distance behavior, which has 
been conjectured in various contexts to be universal. Indeed we find 
that it is totally independent of both Hq and P{'^)- 



^ Unite propre du Centre National de la Recherche Scientifique, associee a I'Ecole 
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1 Introduction 



We consider a Hamiltonian H — Hq + (p with a deterministic Hq and a 

random (/?, in a finite dimensional space in which Hq and (p are Hermitian 
NxN matrices; we are interested in the large N limit. The average number of 
eigenvalues in an interval is known to depend sensitively upon the spectrum of 
Hq and upon the probability distribution P{'p) of the disordered potentiaP^. 
The one-point Green's function and the two-point connected Green's function 
are given by 

«w =< > 

11 11 

- < TrTr 77 >< i^Tr > (1.2) 

N z-H N w-H ^ ' 

We obtain the correlation function between two eigenvalues from G2c{ZiW) 



1 

G^ci}^ + ie, /J, — ie) — G2c{^ — ie, 11 + ie)] (1.3) 



P2c(A,Ai) = - — [G2c{>^ + ie,n + ie) + G2c{\-ie,n-ie 



In the simpler case in which Hq is zero, it was proved in earlier papers^^'^^ 
that the resulting expression for the connected two-point Green's function, 
is independent, up to a scale, of the probability distribution of the random 
potential ip for arbitrary values of A and /i. One word of clarification is 
needed at this stage. While in ref[2] we computed p2c(A, yu) directly, in ref[3] 
we computed G2c{zi,Z2) by letting N go to infinity first, with complex zi 
and Z2- To obtain the correlation function we then let the imaginary parts 
of Zi and Z2 go to zero. By this procedure we obtain a correlation function 
P2c(A, //.) in which a smoothing of the eigenvalues in a range much larger than 
1/N, but small compared to the average point spacing, has been performed. 
This smoothing eliminates all the non-universal oscillations which are present 
in p2c(A,/i) if we had let the imaginary parts go to zero first^^ before letting 
N go to infinity. 

For given non-zero Hq, we could first question whether or not p2c depends 
upon the probability distribution of ip and next, how the result depends upon 
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the spectrum of Hq. We shall answer these questions in several steps: first 
return for definiteness to the pure random case in which Hq = 0, which 
was studied earlier^). Then we consider Hq non-zero and a Gaussian P{^), 
on the basis of earlier work on the subject.^-* We consider next the non- 
Gaussian case with arbitrary Hq and derive first a closed expression for the 
one point Green's function. We then apply this closed expression for two 
specific examples in which the eigenvalues of Hq are restricted to ±1 and in 
which Hq has uniformly spaced levels between —1 and +1. Further we derive 
the two-point connected Green's function G2c{z,w). We find that the short 
distance behaviour is universal. 

The universality of the short distance correlation is perhaps expected 
from the following heuristic argument. Imagine switching on adiabatically 
the potential (p. Starting from some initial eigenstate of Hq, the potential 
will induce a succession of transitions between those eigenstates. For finite 
time it is expected that this process depends sensitively on the nature of 
the spectrum of Hq. However for large times the system will explore all the 
eigenstates of Hq and one might wonder, or may be even suspect, whether 
the limit is insensitive to Hq. Therefore it is natural to consider the limit in 
which the eigenvalues A and /i are near each other, and we prove that in this 
limit p2c is universal. 



2 The simple Wigner case 

We return to the simplest case in which Hq is zero and P{(p) is a Gaussian 
distribution characterized by 

< ipij > = 

< (Pij(Pkl > = J^^jk^il (2-1) 

Much is known about this Gaussian case of course^^'^^ . We have found that 
the correlation function in the large N limit^^ may be written in the compact 
form 

N^G2c{z,, Z2) = -^^Log[l - G{z,)G{z2)\ (2.2) 
It is convenient for later generalization to define 

u{z) = z- G{z) (2.3) 
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given explicitly in this simple case by 




(2.4) 



Since G{z) = l/[z — G{z)], ( |2.2| ) is expressed also by 




(2.5) 



It is useful to introduce also the polarization function 



Uo{zi,Z2)=Giz,)G{z2) 



(2.6) 



We are working here in units in which v, the width of the Gaussian distribu- 
tion, is one; it is easy by dimensional analysis to set it back to any desired 
value. 

There are several ways of deriving ( |2.2| ). In ref[2] it was obtained through 
a study of the asymptotic form of the relevant orthogonal polynomials of high 
order. In ref[4] it was shown that it can be recovered simply by summing 
the diagrams obtained by expanding in powers of v: they consist of planar 
connected diagrams made of ladders and crossed ladders, rainbow-like vertex 
corrections and full planar propagators. 

The form that we give here is not the one which appears in ref[2], but 
it is easy to carry out the differentiations and verify that it is same. By 
differentiating (|2^ ), we have 




(2.7) 




P2c(A,/i) 



1 1 A-Xjj 



(2.8) 



2NV (A - /i)2 [(4 - A2)(4 - /i2)]i/2 



For A close to n it has a singular behaviour 



P2c(A,/i) 



1 



(2.9) 
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It is instructive to compare it with the exact expression^-', i.e. without 
smoothing. In the limit X — n goes to zero with 

X = 27rN{X- fi)p{hx + fi)) (2.10) 

fixed. In this limit the exact correlation function becomes^-* 

^-^^'^)"-^(X^ (2.11) 

The smoothing has thus replaced sin^x by 1/2, thereby producing a double 
pole. 



3 General probability distribution without Hq 

We now go to an arbitrary probability distribution for the disorder. The 
density of eigenvalues is no longer semi-circular; we limit ourselves here to 
spectra which extend over a single segment of length 4a on the real axis. We 
choose the origin at the center of this segment. This is automatic when the 
probability distribution P(v5) is even; otherwise one would have to shift the 
eigenvalues. Then we know from earlier work^-'''^^'^'' that, rather remarkably, 
TSf) is still valid with 

u{z) = -[z + Vz^-Aa^] (3.1) 



2' 



and hence 



P2c[X, fl) = -2iV2^2 _ [(4^2 _ ^2)(4^2 _ ^2)]l/2 ^'^■^) 

The non-Gaussian terms of P{^p) have stretched the support of the spectral 
density from [-2, 2] to [-2a, 2a], and also deformed the density of eigenvalues 
from Wigner's semi-circle to a polynomial multiplying a square-root^^. How- 
ever, the only effect of these non-Gaussian terms on G2c(-2i, -^2) and p2c(A, p) 
is a simple rescaling 2:1^2 — ^ ^ ~^ a' 1^ followed by the multiplying 
G2C and p2c by {l/af- 
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To derive these results we had to proceed in two steps. First^-*, we con- 
sidered a probabihty distribution 

P(y,) = lexp[-iVW(v9)] (3.3) 

in which V{ip) was an arbitrary polynomial of ip. In these cases we could 
diagonalize integrate out the unitary group, and handle the eigenvalues 
of (p by orthogonal polynomial techniques^-*'^"^ (or Dyson's 2D electrostatic 
approach^)'^^^). It was then extended to non-unitary invariant measures, or 
equivalently to a measure in which TrV{(f) was replaced by a sum of arbitrary 
products of traces of powers of ip^\ 

If we were to try again to calculate by expanding the non- Gaussian 
terms of the measure and summing the planar diagrams, we would have to 
deal with expectation values of the type < Tr(y9"Try)'" > expanded in powers 
of the non-Gaussian terms. An equivalent way of stating the results above 
is the following: Call a„,m the Gaussian connected expectation value 

an,n. =< Tnp'^TTip"' >c,o . (3.4) 

For arbitrary n and m it is rather cumbersome to determine these numbers. 
They are of order one in the large N limit, whereas the disconnected part is 
of order A^^. If needed explicitly, we could obtain them by expanding their 
generating function 

N'G^M, Z2) = E «n,^^r"-'^2~'""' (3.5) 

given by ( p.5| ), in inverse powers of the z's. 

We now consider the same expectation value < Tryj^Try)™ >c for the non- 
Gaussian case. The result of (3.2) is equivalent to the following statement 

< Try.^Tr^'" >,= a"+"^«„,^ (3.6) 

in which 2a is again the end point of the support of the spectral density. 

A direct derivation of these identities is non-trivial. For specific values 
of m and n they may be extracted from what is known about G{z). See 
Appendix A. 
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4 One-point Green's function with non-zero 
Ho 

We now have a given unperturbed Hamiltonian Hq, with eigenvalues €{, i = 
1,...,N, with large as always. The one-point Green's function, or the 
average resolvent, was first found by Pastur^^^ for the Gaussian distribution 
P{ip)- Define the deterministic unperturbed one-point Green's function 

Goiz) = ^Tr 







^E— (4.1) 



Then the full one-point Green's function 



(4.2) 



is obtained by solving the implicit equation 

G{z) = Go{z - G{z)) (4.3) 
In other words, the one-point Green's function is determined by 

= < ^Tr -J > 

N^z-ei-G{z) ^^-^^ 

In the simplest case, in which Hq = 0, the above equation becomes a sim- 
ple quadratic equation for G, from which one recovers immediately Wigner's 
semi-circle law for the imaginary part of G, which is proportional to the den- 
sity of eigenvalues. It is easy to give a diagrammatic proof^-* of ( [4.4|) . In 
the large limit, the planar diagrams are the rainbow diagrams whose sum 
leads directly to (gj). 

The self-energy part of the one-point Green's function G{z) is simply 
given for a Gaussian distribution P{(p) by 

E{z) = G{z) (4.5) 
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For the non-Gaussian case, the self-energy part T^{z), defined by 

is of course no longer simply equal to G{z). 

Here we will restrict ourselves to the non-Gaussian probability distribu- 
tion for 

P(^) = i exp[-iV(lTr¥p2 + <^Tr^^)] (4.7) 

To treat the non-Gaussian term gTrip^, it is convenient to consider an "equa- 
tion of motion" obtained by shifting the random matrix by an arbitrary 
matrix (p to (f + e in the integral which gives < {1/ {z — Hq — >■ The 
invariance of this integral under this shift tells us that to first order in e, the 
coefficient of must vanish: 

. 1 . , 1 . 
^^z-Ho-ip'^'^^z-Ho-ip'"''^ 

N < ( ^ )^JiVmn + 4^(<^^)mn) >= (4.8) 

Z — Ho — If 

Setting n = i and m = j, and summing over these indices, we obtain the 
following equation with H = Hq + ip, 

By the definition of two-point connected Green's function, we write above 
equation as 

NG\z) + NG2c{z, z) T:(^-(^) > +4^ < > (4.10) 

Z — 12 Z — 12 

where G2c{z, z) is order of 

The first term on the right hand side of this equation may be rearranged 
by writing ip — [ip -\- Hq — z) -\- {z — Hq) : 

^<Tr{ ^ ip)> = -1 + 1 < IV i {z-Ho)> 

N z — Hq — ip N z — Hq — (f 

N^z-ei- E{z) 
= J^{z)G{z) (4.11) 
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with the self-energy T,{z) defined above. Repeating the same procedure, we 
write the second term on the right hand side of (4.6) 

<IV( i ip')> = 

z- Ho- if \ 

+ < Tr[ ^ J z - Ho) ip{z - Ho)] > (4.12) 
z - Ho - (p 

We have thrown away the term < TYip{z — Hq) > since it is zero by being 
odd in (f. 

We will now consider the relevant Feynman diagrams. As in ref[4] we 
find it convenient to use the language of large N QCD to describe the di- 
agrams. Let us introduce the the gluon-quark-quark vertex function T{z) 
defined to include gluon self-energy corrections. In other hands, it is one- 
particle-irreducible in the quark lines, but not in the gluon line. Then we 
can write the second term in the equation above as 

Y112i^ - ^i){^ - ej) < (- 77 > 

i J z-Ho-if 

- i^^^^'-'^^^'-'\z-e.-m)^-e,-m)^^'^ 

= NV{z){l + i:{z)G{z))^ (4.13) 

As shown in Fig. 1, in the planar limit the vertex function T{z), the 
(quark) self-energy S(2;), and the Green's function or quark propagator G{z), 
are related by 

S(z) = nz)G{z) (4.14) 

Using this fact, and putting everything together, we find that G{z), for arbi- 
trary non-zero Hq and for a non-Gaussian probability distribution for ip, is 
determined by 

G^z) = n^)G{z) - Agx + 4^§||(1 + n^)G{z)f (4.15) 

where we denote < Tr(^^ > by x. This quantity x does not depend on Hq 
and is simply related to the end point a{g) of the spectrum for the Ho — Q 
case by 

^ = y(4-a^) (4.16) 



where a{g) satisfies the following equation^) 



12^a* + = 1 (4.17) 



Thus we have accomplished our goal of finding the equation ( 4.15| ) which 



determines the one-point Green's function G{z) in general. Note that the 
self-energy S which appears in (|4.15| ) is itself related to G{z) as described in 



([4.6|) . These two coupled equations are not explicitly solvable for an arbitrary 
Ho. 

We now consider two simple examples, specified by the eigenval- 

ues ei of Hq: i) ei = e for i = 1, N/2, and = — e for i = {N/2) + 1, N, 
with N even, ii) the are uniformly spaced over an interval from —a to a. 

For the first case i), we have from 



G = ^i J y ) (4-18) 

2 z — e — h z + e — 1^ 



In the case g = 0, ( |4.15| ) tells us simply S = G and thus we obtain the cubic 
equation 

{z-GyG+il-e^)G-z = (4.19) 

The imaginary part of G{z) is represented in Fig. 2 by a solid line for the 
particular case e = 1. The end-points of the spectrum occur at z = ±\/27 /2. 
The point z = becomes also the end point for g = 0. For g 0, we 
consider the following representative values for €=!,&) g = l/2,x = 11/27, 
h) g = —1/48, X = 4/3, and solve numerically. We represent the imaginary 
part of G{z) in Fig. 2 by a dotted line and by a broken line, respectively. 
We have a fifth order polynomial equation for E by substituting ( ^4.18| ) to 
( [4.15|) . At the critical value g = —1/48, we find the imaginary part of S and 



G, vanish at Zc = 3.225 with the exponent 3/2, 

lmG{z) ~ {zc - zfl^ (4.20) 

which shows the same singularity as the if o = case at the edge.^^'^'^-* 
In the second example ii), we have from 



z = J:{z) + a coth{aG{z)) (4.21) 

We consider further the case a = l.The imaginary part of G(z) is evaluated 
numerically, and represented in Fig. 3. At the critical value g = —1/48, the 
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imaginary part of G{z) shows a same singularity at the edge as before. Note 
that the value g = —1/48 always remains critical, independent of Hq. This 
is due to the singular behavior of x given by (4.16) near Qc = —1/48, which 
is independent of Hq. We have considered the non-Gaussian probability 
distribution given by (|4.7| ). It is easy to apply the present method to other 
non-Gaussian P{ip) case. 

5 Two-point correlation function with non- 
zero Hq and (p Gaussian 

We now finally come to the two-point connected Green's function, or cor- 
relation function, G2ciz,w). In an earlier work^-* we have shown that for a 
Gaussian distribution P{(p) we can determine G2c{z,w) as follows. First, we 
define 

For the Gaussian case, the self-energy S = G{z). For notational convenience 
we also define the "circle-product" 

1 ^ 

G-{z)oG-{w) = -Y.g-{z)gT{w) (5.2) 
i=i 

Then in the large N limit, we found 

. G''iz)oGHw) 



N^GoAz^w) 



+ 



^l-G{z)oG{w) 
{G\z) o G{w)){G{z) o G\w)) 
\l - G{z) o G{w)f ^ 
1 1 
l-G{z)oG{z)l-G{w)oG{w) 



(5.3) 



This result was derived by summing all the relevant planar diagrams, consist- 
ing of generalized ladders, with arbitrary cyclic permutations of the rungs, 
with fully dressed planar propagators and planar, i. e. rainbow-like, vertex 
corrections. (We have again set f , the width of the random distribution, to 
one). 



11 



Surprisingly enough, we can show that this rather comphcated expression 
can again be written simply in the form 

iV'G'2e(^i,^2) = -^^Log(l-G(zi)oG(z2)) (5.4) 

Note that this expression has the same form as ( p.2|) but with the ordinary 
product replaced by the circle product. Defining u{z) = z — G{z) as before, 
and using the identity 

7 — in 

l-G{z)oG{w) = —- — (5.5) 

u{z) — u[w) 

we can also write 

AT-in ( \ J ^u{zi)-u{z2) , . 

^ = ^^"^^ z,-z, ^ ^'-'^ 

just as in (p75|). 

To derive these compact representations, we simply differentiate (|5.4|) to 
obtain the previous result ( p. 3| ), using the identity 

(5.7) 



dG{z) G{z) o G{z) 



dz 1 - G{z) o G{z) 



6 Two-point correlation function with non- 
zero Ho and (p non-Gaussian 



The remarkable existence of such a compact representation as ( ^.61) naturally 
prompts us to ask whether this representation also holds for a general non- 
Gaussian distribution P{(p)- In this section, we answer this question partially 
by studying G2c{z, w) for the distribution 

P(^) = i exp[-iV(iTr^2 ^ ^^^^a^^ ^g_^) 

to first order in g. 

As contributions of order g, there are six different classes of diagrams, 
with their typical representatives shown in Fig. 4. The contributions of 
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these six classes of diagrams to G2ciz,w) are 

. G{z)oG{w) 



Da = 



dwdz'l-G{z)oG{wy 



92 ^ G{z)oG{w) , ^2, 



92 . G(z) G2(z)oG(u^) 



+ 



^ a^a^ ' 1 - G{z) o 1 - G(^) o G{w) 
G{w) G^{w)oG{z) 
1 - G{w) o ^(w;) 1 - G{z) o ^(u;)^ 

92 ^ G'3(^) 6-2(^)0 G'(w) 



^1 - G'(z) o G{z) 1 - G'(z) o G{w) 



l-G'(w;)oG'(w)l-(:?(^)oG'(w)J 

To order g(, the one-point Green's function G{z) that appears here can of 
course be taken to be the lowest order Green's function, which we denote by 
Gq{z) (note that this Go{z) is different from the one defined in (4.1), where 
we used it for the deterministic unperturbed one-point Green function.) For 
the sake of notational simplicity, however, we will omit the subscript here. 
In other words, G{z) is determined by 



Summing these terms Da, ...jDf and adding this to the unperturbed term 

dwdz 



N^G2c{z,w) = -^Log(l - G{z) o G{w)) we obtain 



o2 

N^G2ciz,w) = -^^Log{l-G{z)oG{w) 

+ Ag[G{z)oG{w){2 + G^{z) + G\w)) 
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- ''^WoG(,(i^W±^,).0,/) (6.4, 



We note, that to lowest order in g 
G\z)oG{w) 



2/ N „ ^/..A _ G{z) o G{z) - G{z) o G{w) 



w- z- G{w) + G{z) 



G{z) o G{z) - 1 w-z 

+ 77: : ^/ X , (6-5) 



and 



Giz)oGM= , (6.6) 



Then ( p.4[ ) may be written as 

f)'^ y — 7/; 

iV^G,.(.,«;) = -— -Log[( 



X 1 + 4^G(2)G(«;) - 4^ 



9wc?2; z — w ^ G{w) — G{z) 
G{z){2 + G\z)) 



{z-w- G{z) + G{w)){l - G{z) o G{z)) 



. 4 GH(2 + G^H) V 

^(z-«;-G(z)+G(w))(l-GHoG'H);^ ^ ' 

The last two terms of above equation come from parts of diagrams of classes 
(e) and (f). 

We would now like to re-express this result using the exact one-point 
Green's function as far as possible. Using the notation Gq{z) for the one- 
point Green function for non-zero Hq and for the Gaussian distribution ( 
(7 = 0), we find the self-energy T,{z) of the non-Gaussian distribution P{f) 
is given by 

E(.) = G„(.)- «f°W+^ff) +Oto') (6.8) 

1 - Go{z) o Go(^) 

Then the last two terms of (6.7) are absorbed by the expression of the self- 
energy S(w) and T,{z). 

Therefore, we obtain finally the following simple expression, up to order 

9, 

N^G,^{zM = -T^Log[( ' ){l+AgG{z)G{w))] 

fP' 7 — 11) 

-Log[(— — )(1 + ^G{z)G{w))\ + Oig") (6.9) 



dzdw u(z) — u(w) 
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where we now define 

u{z)=z-J:{z) (6.10) 

as the appropriate generalization of (2.3). 

We have thus obtained the same form as in ( ^.6|) except for an extra factor 
of (1 + 4:gG{z)G{w)). This extra factor, however, does not contribute to the 
double pole in p2c{z,w) defined by (1.3). In Appendix B, we show that this 
expression is consistent with the known result for Hq = case. Since this 
factor appears inside log, we separate this term as —4:g{dzG{z)){dwG{w)). 
This term seems to be unconnected part for two-point Green's function. 
Thus if we neglect this factorized term, we have the same form as (2.5). 
The general expression for the connected two point Green's function may be 
written as 

G2c{z,w) = -——Log[l -G{z)oG{w)T^'\z,w)] 

+ ^^-^G{z)G{w)T('\z,w) (6.11) 

The second term represents the diagram class like b in Fig. 4. The fact that 
the terms in ( |6.2|) all collect into the relatively simple form of ( p.9| ) suggests 
to us that some variation of this form may in fact hold to all orders in g, and 
possibly even to arbitrary P{^p)- We do not have a proof of all this tempting 
conjecture at this point. 



7 Discussion 



Thus far, we have considered the correlation function p2c for arbitrary values 
of A and fi, in which case all four terms on the right hand side of (1.3) 
contribute. In the limit in which A — /i tends to zero, one gets a singularity 
which is entirely due to the Green's functions with opposite signs of their 
infinitesimal imaginary parts. Indeed {u{z) — u{w))/{z — w) is not singular 
when w approaches z, with z and w on the same side of the cut. However 
if the imaginary parts of z and w are opposite, and when they approach the 
real axis on the support of the spectral density of the resolvent, the ratio 
blows up at short distance. The residue of the singularity depends upon the 
imaginary part of u, but it disappears when we take the derivative of the 
logarithm. Therefore if A and fi are such that p(A) and p{fi) are non-zero. 
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when ij, is close to A we have a double pole with a universal residue minus 
one: 

N''G2c{\ ±ie,ixT te) = , (7.1) 

(A - fiy 

Consequently, for arbitrary Hq, we have as in the simple Wigner ensemble a 
behaviour 

which, remarkably enough, is totally independent of any specific character- 
istic of the problem. We have proved this for a distribution P{^p) defined 
with a quartic V to lowest order in the quartic coupling, but we are tempted 
to conjecture that this short distance universahty in fact holds to all orders, 
and perhaps even for arbitrary V. 

Wc close with a few concluding remarks. 

(1) Wc have checked that the same result (2.2) holds for a model that 
we have considered recently^^^ , in which random matrices are made of inde- 
pendent random blocks, as when they are attached to a lattice. In that case 
the result comes out immediately from the explicit expression given in that 
work. We obtain 

{CNfG2c{z^,Z2) = -——J2Log{l -ekG{z,)G{z2)) (7.3) 

OZ1OZ2 1^ 

with C the number of lattice sites and the single particle Bloch energies 
as defined in ref[15]. Similarly, in ref[4] we have considered a situation in 
which the random matrix ip depends on a parameter called time. Again, the 
explicit expression given in that work may be written in the form 

N'G2c{zuZ2) = -^£^Log(l-e-"WG(^i)G(z2)) (7.4) 

where u(t) is a function of time defined in ref[4]. In fact the universal form 
discussed in this paper seems to hold in every problem with some random 
features. 

(2) As noted for the simple Wigner case the short distance singularity 
that one finds in (7.1) for the smoothed correlation function is spurious: the 
true p2c(A, n) is finite at short distance, but the smoothing has replaced some 
vanishing numerator by a constant average. This justifies the procedure fol- 
lowed in the literature in computing the universal fiuctuations in mesoscopic 
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systems: the smoothed function correlation is used with some cut-off at short 
distance. An integration by part^^-*, particularly easy since G2C takes the form 
of a second derivative, allows one to return to the integrable logarithmic sin- 
gularity. The short distance cut-off may then be removed. 

(3) Consequently the interesting feature of these results is not so much the 
universal nature of the short distance spurious singularity, but the general 
representation (2.2) of the two-point Green's function. 
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Appendix A: A calculation of < Tup Tnp >c 
Here we will calculate < Tr</7^"Tr(/7^ >c for the probability distribution 
P{<p) = I exp[-7VIV(l(^2 ^ ^^4)] (A.l) 
The one-point Green's function 

satisfies the equation^^^ 

9a2 



G'iz) -iz + 4gz')G{z) + 4gz' + ^-^^^ = (A.3) 



with the end-points ±2a of the spectrum given by 

12c/a^ + - 1 = (A.4) 
Prom G(z) we obtain, after a tedious but simple calculation, 

Then one has 

-|-<4tV'^> = < Tr^'"Tr(^' >c 
og W 



n(n + l)(2-a2)a2"-2^ (A.6) 



n\{n + 2)\ ^ ' dg 

Taking da^/dg — — 12a^/(2 — a^) from (3.10), we obtain as announced 

< Tr(^2n^(^4 a2n+4^^M: = ^^"+4 ^ ^^(^^"TV/ >cO (A.7) 

n -I- 2 (n!)"' 

The method extends to higher powers of m, but becomes very tedious. The 
indirect derivations of ref.[7, 10] are of course easier. 
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Appendix B: Two-point Green's function 



The result of two-point correlation function (|5.4| ) with non-zero Hq and 
with a Gaussian distribution P{(p), is also derived by considering an "equa- 
tion of motion" , which is obtained by shifting the arbitrary random matrix 
(f to ip + e. The propagator {l/{z — Hq — ip)ij{l/{w — Hq — f))ki has the 
corrections of order e by this shift, and they become 



1 . ,1^,1 



( rj")*j( 77") kli,(^nm^mn ~^ ^d^nmi,'^ )mn) (-^•-^) 

Z — II W — il 

By setting n = i, m = j and k = I, and summing over these indices, we 
obtain the following equation, 

< (— Tr ) (— Tr ) > + < — Tr ( )^ > 

- < ^Tr(^-(^ + 4<7y,3)) 1 Tr^- >= (B.2) 

IS z — H M w — H 

The second term of above equation is written by 

<Ti{^—{^^f)> = -A<Tr(^ ^)> 

z ~ H w — H dw z — H w — H 

= -A(gM-gM) (B.3) 

dw w — z 

By the factorization in the large N limit, the first term of (B.2) is given by 

^^Jl^'vhl^'^J^^'^;^^'' = G\^)G{w) + 2G{z)G,^{z,w) 

+ GU^,z)G{w) + Oi^^) (B.4) 

From ( |4.10D and (B.2), we have 



r.^/ / N 9 .G(z) - G(w) 
2Giz)G2ciz,w) + —{- 



dw z — w 

< Tr(^-(yp + Agcp^))TT^— >,= (B.5) 
z — H w — H 
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We consider the two point Green's function G2c{z, w) based on this equa- 
tion. First, we consider the simplest case: Hq = and g = 0. In this case, 
we write the equation for G2c{z, w), 

(2GW-.)G.(.,.) = -A(£(£i^) (B.6) 

Since G{z) — l/{z — G{z)) in the case Hq — and g — 0, we have 

1 G{z) 1 



2G{z) - z G^{z) - 1 v^i^ 
Taking the derivative of w in (B.6), we have 

G2c(z,w) = -^^^—^(G{z)-G{w) + {z-w): ^'^""^ 



(B.7) 



' K./""^ .-1) (B.8) 



2{z - w)2 ' ^/i2T:4 - 4 

This coincides with the previous expression given by (2.7). 
For Ho and g — case, we have from (B.2), 



/r^^/ \ / \ 9 ,G(z) — G(w). 

We diagonahze the matrix Hq with the eigenvalues in the last term. The 
diagrams of this term are classified into four different types (Fig.5). We note 
that if we replace by (e^ + Tj{z) — z) + {z — S(z)) in the all four different 
type diagrams in Fig. 5, we have only an identity without any result for 
G2c{z,w). Namely, in this case we have after calculations, , 

< Tr( i Ho)Ti ^ >e= {z - 2G{z))G2c{z,w) 

Z — Hq — (p W — Hq — (p 
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Inserting this into (B.9), we find that G2c{z,w) cancels out, and obtain the 
following identity equation, 

dw z — w 1 — G{z) o G{w) ow 

To get the expression for G2c{z-,w), we write separately for the diagrams 
of Fig. 5, (a) and (c) as Cj = (cj + — w) + {w — and for (b) and 

(d), ti = [e-i + T,{z) — z) + {z — ^{z)). Then we have the following expression 
by summing the contributions of the diagrams (a), (b), (c) and (d), 

< Tr{ i Ho)Tr 1 {w - E{w) - G{z))G2c{z,w) 

Z — Hq — if W — Hq — (fi 

+ {z-w- E(z) + J:{w))G{z) o G{z)G2c{z, w) 

G\z)oG{w) 



(1 - G{z) o GH)(1 - G{w) o G{w)) 

{G\w) o G{z)){G{z) o G{z)) 
(1 - G{z) o GH)2(1 - G{w) o G{w)) 



(B.12) 



Using the previous identity of (B.ll), we write the equation for G2c{z,w) by 
(B.IO) and (B.12), 

{G{z) - G{w) -z + w){l- G{z) o G{z))G2c{z, w) 
1 



(1 - G{z) o - G{w) o G{w)) 

X [G{z) o G\w){l - G{z) o G{z)) - G^{z) o G{w){l - G{z) o G{w))] 

(B.13) 

It is easy to sec that this equation leads to the previous result (5.3), which 
was derived by a diagrammatic analysis*^^. To see this equivalence, we write 
the following two quantities as 

G\z)oG\w) = i-E ' ' 



N^{z-e,- G{z)y {w-ei- G{w)y 
G'^{z) o G{w) - G{z) o G'^iw) 



w - z- G{w) + G{z) 



(B.14) 



(7-(.)o(;(.) = ^(-)°^(-)-^/-);fH (B.15) 
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and replace these quantities in (B.13), then we obtain the same result as 
(5.3). 

For the non-Gaussian distribution P{(p), we obtain the expression of 
G2c{z,w), up to order g, as (6.9). When we put Hq = 0, this expression 
should be consistent with the known result, discussed in section 3. We will 
show in the following that indeed the expression (6.9) becomes consistent 
with the known result when Hq = 0. 

When i^o = 0, we have by the scaling z ^ z/a and w — > w/a, from the 
expression of G2c{z,w) for the Gaussian distribution P{(fi), 

ffi 7 111 

G2c{z,w) = - Log(l - Go(-)Go(-)) 
ozow a a 

z — w 

Now we put i^o = in the expression of (6.9). Then all circle products 
become usual products, and we have also the following equation, 

G(.)GH = ' ^'"l , , (B.17) 

^ ' ^ ' w- z- S(w) + T.{z) ^ ' 

Then we write (6.9) in the following form, 

2; — w 

G2c{z,w) = - Log' 



dzdw °^z-w + E{w) -J:{z)+ 4:gG{z) - 4c/G(w) ' 

(B.18) 

where the self-energy T,{z) is given as (6.8), 

^ r< ^ \ 8gGo{z) 4:gGl{z) /-r 1n^ 

with 



Goiz) = (B.20) 

We expand the quantity a in (B.16) as = 1 — 12g + 0{g'^), and compare 
(B.16) and (B.18). Then we find that (6.9) is consistent with (B.16) when 
we put Hq — 0. 
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Figure caption 



Fig.l The self-energy T,{z) is expressed by the product of the vertex part 
T{z) and G{z). 

Fig. 2 The imaginary part of the one- point Green's function G{z), for g = 
0, the solid line. The dotted line represents the imaginary G{z) for 
g = 1/2, and the broken line represents for g — —1/48. In all three 
cases, the eigenvalues of Hq are ±1. There is a symmetric counterpart 
for z <0. 

Fig.3 The imaginary part of the one-point Green's function G{z), for the 
case in which the eigenvalues of Hq are given uniformly between — 1 

and 1. The solid line shows the case of 5^ = 0, and the dotted line 
corresponds to g = 1/2, x = 11/27 and the broken line shows the case 
^ = -1/48, a; = 4/3. 

Fig. 4 The diagrams of order g for G2c{z, w). The broken line represents the 
ladder of ip propagators. 

Fig. 5 Four different diagrams of < T^^i z-Ho-v ^o)^ w-Ho-v dotted 
line means the eigenvalue of Hq. 
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